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The topic treated is that of finding a reproducible, plausible and computationally simple method 
of selecting a discrete frequency distribution with a prescribed ranking of its components and pre- 
scribed upper and lower bounds on these components. The problem is shown to be tractable when a 
minimax error selection criterion is employed, and "error" is measured by maximum absolute devia- 
tion among components. In this case one obtains a linear program of a special form admitting explicit 
solution. 

Key words: Linear programs; mathematical models; minimax estimation; operations research; 
probability distribution. 

1. Introduction 

A discrete univariate probability distribution will be represented here as a probability n-vector, 
i.e., an /i-vector x whose components x\ are nonnegative and sum to unity. Let P be some class 
of such distributions, describable as a closed subset of x-space. Then one can pose the problem 
of minimax error selection of a member of P, i.e. of choosing xeP to minimize 

F(x) = max{</(x, y):yeP} (1.1) 

where d is some metric on ra-space. The particular metric 

c/(x,y) = msiXi\xi — y t \ (1.2) 

will be employed, essentially because of its tractability for what follows. 

In [l], 1 this selection problem was solved for the case in which the class P of admissable dis- 
tributions was specified by componentwise bounds, i.e., 

O^L^x^U (1.3) 

where L and U are given n-vectors. (It was also solved for a more general "weighted" version of 
the metric d.) In [2] it was solved for the case in which P was described by a given componentwise 
ranking, 

xi^x 2 ^ . . . ^x„. (1.4) 

The present note gives a solution method for the case in which both types of constraint are imposed, 
i.e., P is specified by (1.3) and (1.4) together. 

Some preliminaries are presented in section 2. The formulation of the solution method is then 
begun in section 3, and completed in section 4. 



1 Figures in brackets indicate the literature references at the end of the paper. 
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2. Preliminaries 

This section considers three subproblems which will arise later. 

PROBLEM I: Given n-vectors A and B, and number S, find an ra-vector y such that 

A^y^B, . (2.1) 

%m = S, (2.2) 

y, ^ y 2 ^ . . . ^ y n . (2.3) 

This problem is treated in [3] (see "Problem IV" in section 2 of that paper), where a solution 
method is given, and the following two conditions are shown to be jointly necessary and sufficient 
for the existence of a solution. First, 

Ai^Bj for i^j. (2.4) 

In terms of the nondecreasing sequences {^}i l and {i?*}? defined by 

A[ — maxj<:j/4j, B[ = minj&iBj, (2.5) 

the second condition is 

ZiA'^S^ZiBi (2.6) 

PROBLEM II: Given ra-vector Z and number S, find the minimum value z* of z such that z ^ and 

2*max{0, Zi~z}^S. (2.7) 

This problem is a special case of one treated in [3] (see "Problem II" in section 2 of that paper), 
where a solution method is given, and 

S ^ (2.8) 

is shown to be a necessary and sufficient condition for a solution to exist. 

Problem III: Given ^-vector Z and number S, find the minimum value z** of z such that 
z ^ and 

S/min {Z/,z} ^S. (2.9) 

This problem too is a special case of one treated in [3] (see "Problem III" in section 2 of that 
paper), where a solution method is given, and 

ItZi^S (2.10) 

is shown to be a necessary and sufficient condition for a solution to exist. 

3. Analysis 

We turn now to the original problem, concerning the constraint set 

P= {x:0 *s L ^ x *S U, 2fXi=l, (3.1) 

Xi ^X 2 ^ • . • ^*n}. 

Recall that the problem is to find xeP to minimize 

F(x) = max {maxi\xi — yt\:yeP} . (3.2) 

As a first application of the material in the preceding section, consider the question of whether 
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the constraints are consistent, i.e., whether P is nonempty. The existence of a yeP is an instance 
of Problem I, with 

A=L,B = U,S=1. 

Thus the necessary and sufficient conditions for consistency are that 

Li^Uj fori^j, (3.3) 

and, in terms of the quantities 

Li = maxj^/Lj, U\ = min,>j Uj , (3.4) 

that 

2,tH*l*2,iUi. (3.5) 

These conditions will be assumed satisfied from here on. Note that P can be rewritten 

P={x:0^L'^x^U\ £i*i=l, (3.6) 

Xi ^ Xi ^ . . . ^ x n } 

where L' and U' are the vectors with respective components L/ and [//. 

Next, set 

Mt = maix {jt/ixcP}, (3.7) 

Mr = mm {xiixeP}; (3.8) 

the determination of these quantities will be discussed later. As in [1] and [2], we have 

F(x) = maxy max, max {y,— Xu X\ — y,} 

= max/ max {max y (y,— #,■), max y (jc,— y,)} (3.9) 

= max, max {Mf—xu Xi — Mf}. 

By (3.9), the minimization of F(x) is equivalent to the following linear program; choose number 
z and vector xeP, to minimize z subject to the further conditions 

z^Mf-xu z^Xi-Mr (alii). 
For fixed z, the conditions on x read 

max {L' i ,Mt-z}^x i ^ min {U'^Mj + z} (alii), (3.10) 

2i*i=l, (3.11) 

Xi^X 2 ^ . • • ^X n . (3.12) 

Thus the objective is to find the smallest z for which such an x exists (and to find one). 

Now finding an x to satisfy (3.10) through (3.12) is an instance of Problem I, with S=l and 



^•=max {L;, Mf-z}, £,= min {U\, Mj- + z}. 



Since sequences {/>,'}'/ and {£/!}? are nondecreasing, while definitions (3.7-8) and the definition 
of P imply that {Mf}j and {Mj} 7 J are nondecreasing, it follows that {Ai}\ and {Bi}'[ are non- 
decreasing. Hence in this case, the conditions for the existence of a solution to Problem I become 

max {L^Mf-z}^ min {£/;, Mj + z) (all i), (3.13) 
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2/ max {/>;, Mt-z}^l, (3.14) 

Si min {{J I, Mf + z} ^ 1. (3.15) 

We seek z min , the smallest value of z satisfying these three conditions. The redundant condition 
z ^ can also be imposed. 
Let us set 

z°= least z satisfying (3.13), 
z* = least nonnegative z obeying (3.14), 
z** = least nonnegative z obeying (3.15). 

Then we have 

z min = max {z°, z*, z**}; (3.16) 

this is the case since the left-hand and right-hand sides in (3.13) are respectively nonincreasing and 
nondecreasing functions of z, while the sums in (3.14) and (3.15) are respectively nonincreasing 
and nondecreasing. 

Since L\ ^ £/■, the determination of z° involves the conditions 

L\ ^My + z,Ml-z^ U'i, Mf-z^Mi + z, 

and leads to 

z () =max; max {Li-Mr,Mt-U' i9 (Mf-Mf)/2}. (3.17) 

Next, (3.14) is equivalent to 

Si max {0, (Mt-Ll)-z}^l- 2/L| , 

so that determining z* is an instance of Problem II with 

Zi = Mt-L'i, S=1-2;L;. 
Finally, (3.15) is equivalent to 

Si min {Ul- Mr, z} ^1-I/M]-, 
so that determining z** is an instance of Problem III, with 

Zi = Ui-Mi, S = l-liMi. 

4. Determination of Mf and Mf 

It only remains to discuss the determination of the quantities M + and Mj defined by (3.7) and 
(3.8). For this purpose, consider the conditions under which 



These conditions are 



x=(yi, . . .,y,-_i, xu y u • . .,r w -i)eP. 

L'i^Xi^U'u (4.1) 

L'j ^ y } ^ min {C/j, xt} for ; < i, (4.2) 

max {L' }+x , x«} ^ yj ^ U' j+i for ; 2* i, (4.3) 

2jyj=l— *«, (4.4) 
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Yi^Ti^ . . . ^Jn-i- (4.5) 

Now the determination of an (n — l)-vector y, satisfying (4.2-5). is an instance of Problem I 
with n — 1 replacing //, with S= 1 —xu and with 

Aj = L'j and Bj = min {c7j, r,} for / < /. 

Aj=xn3L\ {/>j f |. .v,} and Bj=Uj +l for j ^ /. 

Since {L/}," and {c//}," are nondecreasing, and (4.1) holds, the sequences {^j}"" 1 and {/ij},"" 1 are 
also nondecreasing. Thus the conditions in this case, for Problem I to have a solution y, are 

Lj ^ min {Uj , xi} for j < i, (4.6) 

max {£}+!,*,•} ^ f/j +1 fory^i, (4.7) 

2j<iLj + 2j>< max {Lj, x;} ^ 1— ^, (4.8) 

2j<i min {Uj,Xi} + 2j>iUj& 1-Xi. (4.9) 

Conditions (4.6) and (4.7) are automatically satisfied. Condition (4.8) can be rewritten 

Xi + Xj>i max {Lj,Xi} ^ l-Xj<*L}; (4.10) 

the left-hand side is a continuous increasing function of x/, bounded neither above nor below, 
which therefore equals the right-hand side for a unique value xf of x\. Condition (4.9) can be 
rewritten 

*f + 2j<« min {£/],*,} *l-2j>iC/j; (4.11) 

again the left-hand side is an increasing function of*/, which equals the right-hand side for a unique 
value xj of Xu It follows that (4.8) and (4.9) are equivalent to 

Xi ^ Xj ^ x] . 

By combination with (4.1), this yields 

Mt = min{Ul,xt}, (4.12) 

Mr = max {L,\ xy}. (4.13) 

Now the determination of xi and xj must be discussed. 
To determine xf , form the quantities 

LZ=(k-i)Ll. + l Jsk .L'j (k>i). 
The sequence {£*}? is nondecreasing, since 

L^-LZ={k+\-i){L' k+l -L' k )^0. 
We have 

and so there is a last /re{i, i+ 1, . . ., n} for which 

For this A, it is readily verified that 

xt=[i-2j<iL}-l,j >k L}y(k+l~i) ( 4.14) 
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satisfies xfe[L L '■ •> ^a-+i ] • • • or x t ^ L' n if k— n . . . and also satisfies (4.10) as an equation. 
To determine xj , form the quantities 

Ut=(i-k)U' k +3, j «kU' j (k^i). 

The sequence {[/* }/ is nondecreasing, since 

Ut-Ut-i = (i-k + l) {U k -U k _, ) ^ o. 

We have 

i-2,>, u}<i.i«iu;=ut, 

and so there is a least /ce{l, 2, . . . , i} for which 

1-2^ ffj<0J. 

For this A:, it is readily verified that 

* 7 =[l-2 j>( £/;-2j<* l/jl/(£+i-4) (4.15) 

satisfies *re [£//._!, W] • • • or %,~ ^ U[ if fc=l . . . and also satisfies (4.11) as an equation. 

5. References 

[1] Goldman, A. J., and Meyers, P. R., Minimax error selection of a discrete univariate distribution with prescribed com- 
ponentwise bounds, J. Res. NBS 72B (Math. Sci.), No. 4, 263-271 (1968). 

[2] Goldman, A. J., Minimax error selection of a discrete univariate distribution with prescribed componentwise ranking, 
J. Res. NBS 72B (Math. Sci.), No. 4, 273-277 (1968). 

[3] Goldman, A. J., Minimax adjustment of a univariate distribution to satisfy componentwise bounds and/or ranking, 
J. Res. Nat. Bur. Stand. (U.S.), 73B (3) 231-239 (1969). 

(Paper 73B3-302) 



230 



